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We describe a class of exactly-solvable models of interacting bosons based on the algebra SO(3,2). 
Each copy of the algebra represents a system of neutron and proton bosons in a given bosonic level 
interacting via a pairing interaction. The model that includes s and d bosons is a specific realization 
of the IBM2, restricted to the transition regime between vibrational and gamma-soft nuclei. By 
including additional copies of the algebra, we can generate proton-neutron boson models involving 
other boson degrees of freedom, while still maintaining exact solvability. In each of these models, we 
can study not only the states of maximal symmetry, but also those of mixed symmetry, albeit still in 
the vibrational to gamma-soft transition regime. Furthermore, in each of these models we can study 
some features of F-spin symmetry breaking. We report systematic calculations as a function of the 
pairing strength for models based on s, d, and g bosons and on s, d, and / bosons. The formalism 
of exactly-solvable models based on the SO(3,2) algebra is not limited to systems of proton and 
neutron bosons, however, but can also be applied to other scenarios that involve two species of 
interacting bosons. 

PACS numbers: 21.60.Fw, 21.60.Ev, 03.65.Fd 



I. INTRODUCTION 

In the early 1960s, Richardson pj showed how to 
exactly solve the pure pairing model including non- 
degenerate single_j>article orbits for both fermion and 
boson systems p|. In 2001, Dukelsky, Esebbag and 
Schuck jj] showed how to generalize Richardson's solu- 
tion, making use of analogous work by Gaudin [f| for 
spin models, so that the resulting Richardson-Gaudin 
models can now describe the physics of a wide variety 
of strongly-correlated many-body quantum systems gov- 
erned by pairing correlations. Many interesting applica- 
tions of these exactly-solvable models have been reported 
recently, to ultrasmall superconducting grains to nu- 
clear superconductivity , to the nuclear Interacting Bo- 
son Model H, |j| , and to Bose Einstein condensates H3 • 
The methods have been applied both to fermion and bo- 
son systems, invariably yielding useful insight into the 
properties of the complex quantum systems they model 
11]. Furthermore, with slight generalization the methods 
have also proven useful in the description of mixed sys- 
tems involving atoms (either bosonic or fcrmionic) cou- 
pled to molecular dimers in the presence of a Feshbach 
resonance |l2j . 

All of the above models were based on the treatment of 
pairing correlations involving identical particles. The rel- 
evant algebraic structure for describing identical-particle 
pairing is SU(2) for fermions or SU(1,1) for bosons. Both 
are rank-one algebras with just one Cartan generator. 
There has been several efforts to derive exactly-solvable 
Richardson-Gaudin models based on larger rank algebras 
[l3L IT^]. In particular, the rank-two SO(5) algebra de- 
scribing isovector proton-neutron pairing was treated in 
(l5j . More recently, the complete solution of the SO (5) 
Richardson-Gaudin model was presented in ref. |l6l ] . This 



made possible an exact treatment of isovector proton- 
neutron pairing in the presence of non-degenerate single- 
particle levels, which was used in a description of the 
N = Z nucleus 64 Ge in an extended shell-model space. 

In this paper, we describe the first example of an 
exactly-solvable Richardson-Gaudin model of interacting 
bosons that is based on a rank-two algebra. The model 
has the algebraic structure SO(3,2) and is of relevance 
to two-component boson systems. One possible realiza- 
tion is a mixture of 97 Rb atoms in the hyperfme states 
|F = l,M f = 1), \F = l,Mf = -1) H3. While these 
systems have been studied [l|| in mean-field approxima- 
tion, more elaborate approximate methods or the use of 
exact solutions are needed to study effects beyond mean 
field. Here we will concentrate on a different realization, 
namely bosonic models of nuclei with distinct proton and 
neutron degrees of freedom. When only s and d bosons 
are included, the resulting models are restricted versions 
of the proton-neutron Interacting Boson Model (IBM2) 
|l9j . in the transitional regime between gamma- unstable 
and vibrational nuclei. At the same time, however, their 
generality makes it possible to address many issues that 
cannot be discussed within the context of the IBM2. In 
particular, the models can readily accommodate boson 
degrees of freedom other than just s and d, still within 
a proton-neutron framework. As such, they can address 
issues related to the proton-neutron degree of freedom in 
boson models with g bosons and/or / bosons. The mod- 
els can also accommodate F-spin breaking, which can 
thus likewise be studied in the presence of bosons other 
than s and d. 

The outline of the paper is as follows. We begin in 
Section II with a brief overview of the SO(3,2) algebra, 
confirming that it is indeed a rank-two algebra and ob- 
taining its associated integrals of motion. In Section III, 
we show how this formalism can be used to build exactly 
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solvable models of relevance to interacting neutron and 
proton bosons in atomic nuclei. In Section IV we con- 
sider as specific examples models involving sd, sdg and 
sdf degrees of freedom and describe their exact solutions 
in a variety of scenarios. In Section V, we summarize the 
principal conclusions of this work. 



II. SO(3,2) AND ITS CARTAN 
DECOMPOSITION 



In this section, we briefly review the elements of the 
Cartan decomposition of the SO(3,2) Lie algebra and 
show how it can be used to build the set of integrals 
of motion for models based on multiple copies of this al- 
gebra. 

We will be considering a problem involving neutron 
and proton boson degrees of freedom, with creation and 
annihilation operators Zj m and l pm , respectively. These 
operators create/annihilate a boson of angular momen- 
tum /, projection m and type p. 

Now consider the set of bilinear boson operators 
(p, p' — v or 7r) 

P P ~ Z-^, l pm l P m ) 



bin = ■ It , 

bl = (-l)^t ■ It , 

b\ u = {-l) l '%-ll, (1) 

and the corresponding two-boson annihilation operators 
6_n, box, 6+11. Here the scalar product ij, ■ V p , has the 
usual definition 

P P ' P,m b p',-m ■ 



1% 



■D l il ■ I 



(2) 



It can be readily confirmed that [H^Hf] — 0. As 
there are no other generators that simultaneously com- 
mute with both, SO(3,2) is a rank-two algebra and these 
two operators generate the Cartan subalgebra. [Note: 
the factor 1/2 is included in the Cartan operators to en- 
sure that the Cartan- Weyl basis is orthonormal with re- 
spect to the Killing form of ref. [2fT|.] 

The remaining generators are the ladder operators of 
the algebra. and bj are the raising operators and 
correspondingly F7~ and bi^ are the lowering operators. 

Knowing the Cartan decomposition for each I, we can 
obtain the associated set of integrals of motion Ri , follow- 
ing ref. |14| . The complete set of hermitean and mutually 
commuting operators with linear and quadratic terms in 
the rational model can be expressed in the form 



+ £(2 + A)fl? + £ 
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2.9 



X t ■ Xy 

zi> - Zi 



(3) 



where Hf (a = 1,2) are the two Cartan operators of the 
copy I, and 



X l -X v =F r F l , + HfHf,-\ Kl b ^' 



KM 



(4) 

For a problem involving K copies of a rank-two SO(3,2) 
algebra the integrals of motion depend on a total of K + 2 
free parameters, the K coefficients z; and the two param- 
eters, A and g, that enter in the linear term of eq. J3J). 

The exact solutions of the rational model correspond 
to solving the eigenvalue equations 



i?i|*(e Q ,w 7 )) = r/|*(e a ,w 7 )) 
The eigenvalues r; can be written as 



(5) 



For each I, there are four particle-hole operators , three 
pair creation operators, and three pair annihilation oper- 
ators. All told, these ten operators generate an SO(3,2) 
algebra. 

The Cartan decomposition of this algebra is achieved 
by rewriting the set of one-body operators as 



n = ^-(Uiri + u vl + Au nl ) + ^-(2 + A)n ; 



+ 



2 ' zi> — zi 



-^—(fli(U n i> + U vV ) + Slv(U*t + U v i) 



Hi — F z i — ^ |^(4 m ^7rm ^im^m)^' j 



1 



Mi 



M-2 



7=1 



Zl 



a=l 



Zl 



(6) 
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where U p i(p — ir and v) are the number of non-paired p 
bosons (i.e., the seniority) in the level I, and which define 
completely the SO(3,2) irreducible representation of copy 
I. Also, = (21 + l)/2 is proportional to the degeneracy 
of the levels, Mi — N v — J^i U^i is the total number of 
paired n bosons, and M% = M = ^(N^ — J^i U w i + N v — 
J2 t U v i) is the total number of pairs. Finally, ui 1 (7 = 
1, .., Mi), and e a (a = 1, .., M2) are two sets of spectral 
parameters that must satisfy the set of coupled equations 



where 



E 



E 



U,tI — Uvl 
Zl — UJ-f 



zi - e a 



Mi 

E 



Mi 

E 



Cdg — OJ-y 



M 2 

E 



e/3 



e /3 



Mi 

E 



^5 



A 



1 • (7) 



(10) 



c = g X)r>; ^z^i' — (2 + A) e/fi; is a constant, and 1 
is the identity operator. 

The diagonal term in the last line of eq. © depends 

on 

B(U n ,U u ) = 2 U " l +(U^+U u i)Qi-U vl -2U n i 

1 

(11) 

a function of the seniorities, and a constant C = N(4tt + 
6 — N)/4, with = J^i an d N the total number of 
bosons (N = J2iN wl +J2i N »i)- 

The eigenvalues of the above Hamiltonian (2g J^i e i r i + 
c) can be expressed as 



These are generalized Richardson equations appropriate 
to the SO(3,2) algebra. 

Once the free parameters that define the Ri integrals 
of motion have been chosen, any linear combination, 



E^ ' 



(8) 



E e « 



A 



5> 7 -(-F 2 )(-F 2 + l) 



of those integrals of motion also gives rise to an exactly- 
solvable model for the multiple copies of the SO(3,2) al- 
gebra, with eigenvalues obtained from the correspond- 
ing n's. Of course, any constant can be added to the 
Hamiltonian without affecting the exact solvability of the 
model. 

For the scenario under discussion, a system of inter- 
acting neutron and proton bosons in several levels, it is 
obviously desirable that the exactly-solvable models that 
emerge have the appropriate single-particle energies for 
the various boson degrees of freedom. This is accom- 
plished by requiring that the coefficients 77/ that define 
the linear combination of integrals of motion are propor- 
tional to the single-particle energies of the corresponding 
boson levels (77; = 2gei). If we furthermore choose the 
zi parameters that enter in the definition of the integrals 
of motion according to z; = 2e;, the two-body part of 
the model Hamiltonian likewise emerges in a very simple 
form, one in which the interaction matrix elements are 
level- independent. Furthermore, the constant term that 
is added to the Hamiltonian can be chosen so as to leave 
it in a convenient form. 

With these choices, the form of the Hamiltonian that 
emerges for multiple copies of the SO(3,2) algebra is 



H = 2g^eiRi+d 



Nul + AiVjrj 



.l F . F + i( B (U v ,U v ) + C)i, 



(9) 



+Z(B(U W ,U V ) + C) , 



(12) 



where F z = (N^ — N v )/2. The first line is the single parti- 
cle energy of the non-paired bosons. The first term of the 
second line is the energy of the boson pairs, whereas the 
second term is related to F-spin symmetry. It is possi- 
ble to show that in the i^-spin-symmetric limit (A — ► 0), 
both terms in the brackets combine to give: 



lim 
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(-F z )(-F z + 1) 



= -F(F + 1) 



It is especially convenient to rewrite the Hamiltonian 
@ in multipole-multipole form, to facilitate a link to 
the traditional phenomenological Hamiltonians used in 
bosonic descriptions of nuclei. The Hamiltonian @ ex- 
pressed in terms of multipoles can be written as 



H = E ei fa* + ft* + A ^-> 



v+i 

-|E E (-) L (Qtii+QiO<Qfa* + QiiJ > ^ 

Kl' L=\l'-l\ 

where 



Q^ = (iiv P -{-r l -^i'fh 



(14) 
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III. SO(3,2) AND THE PROTON-NEUTRON 
INTERACTING BOSON MODEL 

The multi-copy SO (3,2) formalism described in the 
preceding section has relevance to the proton-neutron 
Interacting Boson Model of atomic nuclei 0]. In this 
model, distinct neutron and proton bosons are intro- 
duced to model collective pairs of identical nucleons. 
The standard version of the model is limited to s and 
d bosons, reflecting the fact that the energetically low- 
est pairs of identical nucleons are invariably those with 
J p = + and J p = 2 + and that they are typically well 
separated from all higher pairs. The resulting model, 
usually called the IBM2, has been used to describe suc- 
cessfully the collective properties of nuclei throughout the 
periodic table. 

Soon after the introduction of the IBM2, it became 
clear that there was a natural connection of the model 
both to the nuclear shell model and to the simpler IBM1 
model with but one kind of boson. 

The relationship to the shell model is contained in the 
association of the bosons with the lowest collective pairs 
of identical nucleons. Thus, it is possible to use input 
from the shell model to help define the Hamiltonian that 
acts in the IBM2 model space. 

The relationship to the IBM1 comes through the intro- 
duction of a quantum number called F spin that distin- 
guishes the two species in IBM2, the neutron and pro- 
ton bosons. In an IBM2 Hamiltonian that conserves 
this quantum number, the lowest states of the system 
are those that are maximally symmetric in F-spin, in 
other words maximally symmetric under the interchange 
of neutron and proton bosons. It is these states that are 
modelled by the simpler IBM1. However, the IBM2 also 
contains states that are not maximally symmetric under 
the interchange of neutron and proton bosons, states that 
are said to have mixed symmetry. In an F-spin conserv- 
ing version of IBM2, these states are decoupled from the 
states of maximal F-spin symmetry. 

With the above as background, we are now in a posi- 
tion to discuss the relationship between the exactly solv- 
able models that derive from the SO(3,2)algebra and the 
IBM2. In particular, if we limit our discussion to two 
copies of the algebra, one involving an I — s boson and 
the other involving an I = 2 d boson, then we arrive pre- 
cisely at a Hamiltonian of the IBM2 form. It is, however, 
not the most general IBM2 Hamiltonian. Rather, as we 
will see shortly it is a Hamiltonian that models the transi- 
tion between vibrational (U(5)) and gamma-soft (SO(6)) 
nuclei within an IBM2 framework. 

The Hamiltonian that derives in a two-copy s-d real- 
ization of SO(3,2) has the form (see eqs. I|13ll4fl . with 
I, I' — s,d only, and taking e s — 0) 

H = e d (fi„ d + N vd + AN^ 



where 

Q 2 p = sld+dl~s p . (16) 

Its eigenvalues can be obtained from eq. i|12[l . 

There are a couple of features of this Hamiltonian wor- 
thy of note. The first is that as mentioned above it 
is not a completely general IBM2 Hamiltonian, as the 
quadrupole operator is restricted to v = 0. Second, 
the Hamiltonian contains a term AN^d, which arises 
if the single-boson splittings between the s and d lev- 
els in the two species (neutrons and protons) are differ- 
ent, and which will in general break F-spin symmetry. 
F-spin symmetry is only preserved if A — 0. Finally, 
the quadrupole-quadrupole interaction that enters is an 
F = tensor, with a definite relation between the in- 
teractions among the different species. There is no ad- 
ditional Majorana interaction in the model, as the term 
gF ■ F in eq. © is absorbed completely into the Jo- 
spin-scalar quadrupole-quadrupole interaction. However, 
when A = 0, i.e. when F-spin symmetry is preserved, 
we can add an additional Majorana term of arbitrary 
strength while maintaining the exact solvability of the 
model. 

Perhaps the key feature of the exactly-solvable proton- 
neutron boson models that derive from the multi-copy 
SO(3,2) algebra is that they can accommodate boson de- 
grees of freedom other than just the s and d. Indeed, 
they can accommodate any number of boson degrees of 
freedom, albeit with the restricted Hamiltonians that can 
be accessed. As an example, the model with s, d and g 
degrees of freedom is of the form (e s = 0) 

H = E ei + ^ + Ai ^) 

l=d,g 

Kl' L=2,...,6 

Note that it contains all multipole interactions from 
L = 2 to L = 6. Once again, all enter as F = ten- 
sors, with no additional Majorana interaction present. 
Here too, however, a Majorana interaction of arbitrary 
strength can be included if A = so that F-spin is con- 
served. 

Likewise, we can accommodate / bosons, either with or 
without g bosons, still in a proton-neutron boson frame- 
work. When considering sd and / bosons, for example, 
the class of Hamiltonians that arise are (e s = 0) 

H = E (l + ^ + A ^ 1 ) 

l=dj 



~{Ql + Ql)-{Ql + Ql) , 



(15) 



-|E E (-) L (QU+QtO<Qh» + QiO • ( 18 ) 

l<l> L=l,...,5 
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1,3, 5 , 7 + 


(0 1 3 0) 
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+ ,3 + ,4+,6+ 
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2+,4 + ,5 + ,6+,8 + 
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0,1,. 


• ■ 5 Fmax 


0+,l+,...,8+ 


(010102) 


0,1,. 


■ ■ ? Fmax 


0+,l + ,...,10+ 



TABLE I: Seniorities considered in the sd and sdg calculations 
and their associated angular momenta, parities and F-spin 
quantum numbers. 



Note that now we only have multipoles from L = 1 
through 5. The L = 1, 3, 4 and 5 multipoles all have 
odd parity; the L — 2 multipole has an even-parity 
component (1,1' = s,d) and an odd-parity component 
(1,1' = dj). 



IV. APPLICATIONS OF EXACTLY-SOLVABLE 
PROTON-NEUTRON INTERACTING BOSON 
MODELS DERIVED FROM THE SO(3,2) 
ALGEBRA 

In this section, we apply the various exactly-solvable 
model Hamiltonians developed in the previous section to 
scenarios of relevance to the collective structure of atomic 
nuclei. 

We first consider a model of s, d and g interacting neu- 
tron and proton bosons, for which the relevant Hamilto- 
nian was given in eq. (|17fl . We compare the results for 
this model with those of the corresponding model with 
just s and d bosons, (see eq. (I15I) ). The analysis is car- 
ried out as a function of the pairing strength <?, for fixed 
single-boson energies, with the results shown in fig. 1. 

The complete Hilbert space of the model can be split 
into invariant subspaces characterized by the seniority 
numbers (U^ijU^i). For each seniority, one has to solve 
a different set of Richardson equations (7J to obtain the 
eigenvalues <|12[) . The angular momenta, parities, and 
F-spin quantum numbers associated with the set of se- 
niorities we consider are shown in table [I] 

The top panel gives the results for the sd model and 
the bottom panel gives the sdg results. All energies are 
with respect to that of the ground state. 

In the case of the sd model, we use td — 1 MeV and 
e s = M eV, and the relevant seniorities are those with 
U g = in table [I] For the sdg model, we again use 
ed = 1 MeV, and e s = MeV, but now include a g 
boson at e g = 1.6 MeV. In this case, we consider all the 
seniorities of tabled The level energies chosen are physi- 
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2 ,4 ,5 ,6 
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FIG. 1: Excitation energies of the sdg Hamiltonian in eq. I|17|l 
(bottom panel) compared to those of the sd Hamiltonian in 
eq. I|15|l (top panel), as a function of the interaction strength g 
and for the seniorities of tableQ] Continuous lines correspond 
to maximal F-spin states, dashed and circled lines refer to 
mixed symmetry states (F = Fmax — 1). Thick and circled 
lines in the bottom panel indicate states with non-paired g- 
bosons (U g ^ 0), for maximal and mixed symmetry states 
respectively. Except for J p — 2 + , only the first excited state 
of each seniority is shown. Only those mixed symmetry states 
up to total seniority two are presented. The single-boson en- 
ergies are given in the text. 



6 



cally reasonable for the various boson degrees of freedom. 
In both calculations, we consider a system of N = 20 
bosons, with N v — N n — 10 (F max = 10). Furthermore, 
we assume A = so that -F-spin symmetry is preserved. 

Note that 5 = corresponds to the precise U(5) limit 
in the sd model and to the corresponding vibrational 
limit in the sdg model. As the pairing strength g in- 
creases, there is a phase transition to a gamma- unstable 
system, occurring at roughly g — 0.08 MeV. 

Some interesting features can be noted in these results. 
Inclusion of the g level in the calculation increases sig- 
nificantly the number of possible J p states (indicated by 
thick and circled lines in the bottom panel of fig. 1, see 
table nj , but otherwise leaves most of the states of the 
sd model relatively unaffected. For a few states, how- 
ever, and in particular the levels denoted + , ms , and 
2\ , there are important and interesting changes that take 
place. As an example, consider the first excited + state. 
In the sd calculation, its energy goes up as g increases; in 
the corresponding sdg calculation its energy flattens out. 
A similar effect can be seen for the ms and 2^~ states, 
which likewise are flattened out in energy. 

To better illustrate this phenomenon and how it arises, 
we present in the bottom panel of fig. 2 the pair energies 
(e Q ) associated with some selected seniority-zero states 
(J p = + ) of the sdg model, as a function of the inter- 
action strength. We also show the associated spectrum 
of levels in the top panel of the figure. 

For the ground state of the system (GS), all pair en- 
ergies are trapped between 2ed and 2e s . Then there are 
a series of levels in which one or more of the pair ener- 
gies are promoted into the energy region from 2e<i to 2e g . 
The corresponding levels are likewise trapped, as their 
energies are simply a sum of the associated pair energies. 
The lowest is a±, with but one pair excited. The next 
is a%, with two pairs excited. And as should be evident 
there are many other such states, in which at least one 
pair is between 2e g and 2ed, but in which none are above 
2e g . For simplicity, however, we only show the lowest two 
levels of this band. 

Next we consider levels in which one pair is excited 
above 2e g . The lowest such level is denoted (5\ in both 
panels. Since the last pair is not trapped, this level goes 
up in energy as g increases. Furthermore, there will be 
a band of levels in which just one pair has been excited 
above 2e g . Once again, however, we only show the two 
lowest states of the band. 

Finally, we show two other bands in the spectrum, de- 
noted 7 and 5. They correspond to two and three pairs 
above 2e g , respectively. They too are unbounded as a 
function of g. Of course there are higher bands as well, 
corresponding to progressively more pairs being excited 
above 2e g . 

As should be clear from the figures and the discussion, 
this unusual behavior for selected levels occurs whenever 
there are more than two boson levels active and a repul- 
sive boson-boson pairing interaction [|J. 

As noted earlier, the trapping of pair energies is not 
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FIG. 2: Excitation energies of the sdg Hamiltonian in eq. I|17|l 
(top panel) for seniority-zero states (J p — + ). The thickness 
of the lines serves as a guide to the eye to distinguish the 
band structure of the spectrum (a, (3, 7,5). Only the first 
two states for each band are plotted. In the bottom panel the 
pair energies (e a ,a — 1, 10) associated with the ground 
state (GS), the first (qi) and second (02) states of the lowest 
excited band, and the first (/3i) state of the second excited 
band are shown (see text for details). All the parameters are 
the same as for the bottom panel of fig. 1. 
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TABLE II: Seniorities considered in the sdf calculations and 
their associated angular momenta, parities and F-spin quan- 
tum numbers. 



restricted to + states. Indeed, the band structure found 
for seniority-zero states is a feature common to all senior- 
ity sectors. However, for the states with seniority differ- 
ent from zero the term (g/2)B(U 1T , U^) in the eigenvalues 
(|12J) also contributes, making the excitation energies for 
such states unbounded as a function of the interaction 
strength g. Nevertheless, these states likewise arise from 
trapped pair energies and thus are also flattened out as 
a function of g. 

For the mixed-symmetry states in the A = limit, the 
term in brackets in l)12ll produces a contribution to the 
excitation energy, 
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FIG. 3: Excitation energies of the sdf Hamiltonian in eq. I|18[l 
as a function of the interaction strength g, for the seniorities 
of table [n] Continuous lines correspond to maximal F-spin 
states, dashed and circled lines refer to mixed symmetry states 
(F = Fmax — 1). Thick and circled lines indicate states with 
non-paired f-bosons ([// 7^ 0), for maximal and mixed F- 
symmetry states respectively. Only the first excited sate of 
each seniority is shown. For mixed symmetry states, only 
states up to total seniority two are plotted. The single-boson 
energies are given in the text. 



thereby explaining why they appear at higher energy 
than the maximal .F-spin states. As was mentioned 
above, however, when A = it is possible to add a Ma- 
jorana term of arbitrary strength, and thus move the lo- 
cation of the mixed-symmetry states at all with respect 
to those of maximal .F-spin. 

We next consider the effect of including an / boson, fo- 
cussing now on an sdf boson model. In fig. 3, we present 
the results for this model as a function of g, for the case of 
e s = MeV, e d = 1 MeV and e f = 1.6 MeV. The senior- 
ities that were considered in this calculation are listed in 
table HH 

Thus, this model is very similar to the sdg model shown 
in fig. 1, with the / boson replacing the g boson at pre- 
cisely the same energy. Here too the second + level 
is trapped, as this is a general phenomenon that oc- 
curs whenever there are more than two boson degrees of 
freedom and a repulsive pairing interaction. Indeed, the 
spectra for the sdf and sdg models as a function of g are 
very similar. The exactly-solvable SO(3,2) Richardson- 
Gaudin models do not couple states with broken pairs, 
which is where differences would have shown up. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have described the first example of 
an exactly-solvable Richardson-Gaudin model based on a 
rank-two algebra for bosons. The example we discussed 
involved multiple copies of the SO(3,2) algebra. We fo- 
cussed on a specific realization involving proton and neu- 
tron bosons subject to a pairing interaction. The mod- 
els that emerged are exactly solvable for non-degenerate 
single-boson energies and for any number of bosonic 
copies. When only s and d bosons are included, the re- 
sulting model is a specific example of the proton-neutron 
Interacting Boson Model (IBM2), limited however to the 
transition from vibrational to gamma-soft motion. Other 
bosonic degrees of freedom, such as the g and/or /, can 
be readily included while still preserving the exact solv- 
ability of the model. Through the use of these models, we 
can address issues related to the role of mixed-symmetry 
states in proton-neutron boson models with many bo- 
son degrees of freedom. We can even study some limited 
features associated with F-spin symmetry breaking for 
these multi-copy boson models. 
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One issue addressed here concerned the role of the g 
boson in proton-neutron boson models. We compared 
the results of sd calculations with those of sdg calcula- 
tions for an F-spin-symmetric Hamiltonian, throughout 
the vibrational to gamma-soft transition region, and saw 
that in the presence of a g boson several low-lying excited 
states have their properties dramatically modified as the 
result of a trapping phenomenon for the pair energies. 

We also studied the role of the / boson, in the con- 
text of an sdf model of interacting proton and neutron 
bosons. Because of some special features of the exactly- 
solvable models that we treat, there were no meaningful 
differences between the sdg and sdf results, except of 
course as regards the angular momentum content of the 
collective states. 

The models that we have developed are not limited, 
however, to proton-neutron boson models of nuclei. Any 
physical problem involving two species of bosons in which 
pairing is dominant can be modelled in this way. In Sec- 



tion I, we noted the possibility of applying these models 
to problems involving a mixture of 97 Rb atoms in the 
hyperfine states \F = l,M f = 1), \F = l,M f = -1). 
The availability of an exactly-solvable model with which 
to study such systems, makes it possible to study effects 
that go beyond the limited mean-field descriptions con- 
sidered to date. 
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